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SUMMARY 

Using  a  North  American  Navion  as  a  flight  test  vehicle,  a  study 
was  made  of  the  practicality  of  experimentally  determining  all  longitu- 
dinal stability  derivatives  uniquely  by  releasing  a  large  weight  in  flight. 
The  solutions  for  the  stability  derivatives  were  attempted  by  matching 
the  transient  response  of  the  aircraft  with  solutions  to  the  equations  of 
motion  generated  on  an  analog  computer. 

Due  to  instrumentation  difficulties  definite  verification  of  the 
practicality  of  this  procedure  was  not  achieved.     However,  the  experi- 
mental results  indicated  that  the  stability  derivatives  could  be  found 
uniquely  by  releasing  a  large  weight  in  flight. 

Unique  longitudinal  stability  derivatives  were  found  from  simu- 
lated data  generated  on  the  analog  computer.     Better  accuracy  for  these 
derivatives  was  achieved  if  the  transient  response  to  an  elevator  step 
was  analysed  in  addition  to  the  weight  release  response. 

It  is  considered  that  this  procedure  for  the  unique  determination 
of  longitudinal  stability  derivatives  is  suitable  for  further  investigation. 
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INTRODUCTION 

In  recent  years,  dynamic  flight  testing  has  received  progres- 
sively increasing  emphasis .     The  objectives  of  dynamic  testing  are 
several  in  number.     However,   receiving  the  greatest  attention  is  the 
determination  of  stability  derivatives  from  the  transient  response  of  an 
aircraft  to  an  applied  forcing  function.     There  are  many  advantages  to 
this  type  of  flight  testing.    Among  the  more  important  advantages  are 
the  short  flight  time  and  relatively  simple  flight  techniques  required. 
In  many  cases,  dynamic  flight  testing  offers  the  only  means  of  finding 
certain  stability  derivatives . 

The  basic  approach  in  using  dynamic  flight  testing  to  determine 
stability  derivatives  is  simply  described.     A  time  history  of  the  air- 
craft's transient  response  to  some  input  forcing  function  is  recorded. 
It  is  assumed  that  the  form  of  a  set  of  linear,   constant  coefficient,  dif- 
ferential equations  describing  the  system  is  known.     The  time  history 
is  then  U6ed  to  evaluate  the  constant  coefficients  in  these  equations. 

A  basic  problem  that  is  present  in  dynamic  flight  testing  is  the 
unique  determination  of  all  longitudinal  stability  derivatives       The  most 
common  procedure  has  been  to  assume  some  theoretical  relation  be- 
tween these  derivatives.     This  report  is  concerned  with  the  unique  de- 
termination of  all  longitudinal  stability  derivatives  without  any  such 


assumptions.      This   is   accomplished  by   using  a  weight  release   as   a 
forcing  function  to  produce  the  aircraft's  transient  response .      This 
relieves  any  linear  dependency  between  the  various  longitudinal  stability 
derivatives. 

This  investigation  was  carried  out  at  the  James   Forrestal 
Research  Center,   Princeton  University,  Princeton,  New  Jersey,  during 
the  1938-1959  academic  year. 


EQUIPMENT 

Airborne   Instrumentation 

The  teBt  aircraft  employed  was  a  North  American  Navion, 
which  i6  a  four  place,  all  metal,  low  wing  aircraft.      The  Navion  is 
powered  by  a  Continental  E-185  engine  rated  at   185  HP  at  2300  RPM 
and  29  inches  manifold  pressure.     A  three  view  drawing  of  the  aircraft 
is   shown  in  Fig.   1,  and  the   general  specifications  are  recorded  in 
Table  I. 

The  aircraft  was  modified  by  the  installation  of  a  Navy  Mark 
51  bomb  rack  which  was  placed  on  the  lower  center  line  of  the  fuselage, 
aft  of  the  nose  gear  wheel  well.     The  bomb  rack  was  installed  so  that 
the  release  mechanism  could  be  actuated  either  electrically  by  a  sole- 
noid or  manually  by  a  cable.    A  566.  5  pound  solid  steel  bar  was  equipped 
to  hang  from  the  bomb  rack.     This  steel  bar  was  9.  5  inches  in  diameter 
and  27  inches  long.    Henceforth  in  this  report,  the  steel  bar  will  be  re- 
ferred to  as  the  bomb.     Figs.  2  and  3  show  the  bomb  installed  on  the 
aircraft. 

Standard  aircraft  instruments  were  used  to  determine  altitude 
and  airspeed.    In  addition,  the  aircraft  was  instrumented  to  measure 
normal  acceleration,  pitch  rate  and  angle  of  attack.     The6e  variables 
were  fed  to  a  telemetering  system  and  recorded  at  a  ground  station. 
Also,  the  instant  of  bomb  release  was  recorded.     Fig.  6  is  a  schematic 
wiring  diagram  for  the  instrumentation  system. 


The  accelerometer  used  was  a  Giannini  Model  25113-1-20  with 
a  range  of    ^1  g,      The  sensitive  axis  of  this  instrument  was   canted 
laterally  at  45°.     With  this  arrangement,  the  accelerometer   will  only 
respond  to  a  45°  component  of  normal  acceleration.     This  effectively 
multiplies  the  range  of  the  instrument  by  the  square  root  of  two  for 
normal  accelerations.     The  full  scale  range  is  therefore    -  1.41   g's. 
An  initial  installation  with  a    -  3  g's  Genesco  accelerometer  was  made, 
but  unsatisfactory  data  was  obtained.     No  other  small  range  accelerome- 
ter was  available  which  could  be  employed  with  the  telemetering  equip- 
ment.    A    i  1  g    accelerometer  canted  laterally  at  45°  will  of  course 
respond  to  any  lateral  accelerations.     However,  the  instrument  provided 
reasonably  good  results.     Full  scale  voltage  for  the  accelerometer  pot  en 
tiometer  was  provided  from  the  28  volt  battery  bus.     The  resistance  of 
the  potentiometer  is  2000  ohms.     A  bias  battery  was  included  to  adjust 
the  steady  state  voltage  to  the  desired  level. 

The  pitch  rate  was  measured  with  a  Honeywell  gyro  Air  Force 
Part  No.   JG  7005  A-24.     This  type  instrument  operates  with  a  rotor 
power  supply  of  single  phase.  400  cycle,    115  volt  alternating  current. 
The  transducer  is  of  the  potentiometer  type  with  a  resistance  of  530 
ohms.     The  axis  of  the  gyro  was  carefully  checked  to  align  with  the  air- 
craft axis.     The  potentiometer  was  supplied  with  five  volts  from  the 
battery  bus.     These  same  five  volts  were  supplied  to  the  telemeter  equip- 
ment as  a  full  6cale  reference. 
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The  angle  of  attack  vane  was  mounted  on  a  four  foot  boom  which 
was  located  on  the  left  wing  tip.    A  balsa  vane  with  a  direct  connection  to 
a  helipotentiometer  was  employed  to  reduce  the  moment  of  inertia  and 
friction  of  the  system.     The  potentiometer  wa6  supplied  with  28  volts  and 
a  bias  battery. 

The  telemeter  unit  was  the  pulse  width  type  and  was  built  by  the 
Applied  Science  Corporation  of  Princeton.     This  unit  contains  43  channels 
and  a  normal  sampling  rate  of  20  per  second  on  each  channel.       Three 
channels  were  used  for  angle  of  attack  and  pitch  rate  while  six  channels 
were  used  for  the  normal  acceleration. 

The  telemetering  unit  is  designed  to  use  five  volts  as  a  full  scale 
reference  voltage.     When  properly  instrumented,  these  same  five  volts 
should  be  the  full  scale  reference  for  each  transducer .     Small  variation  of 
the  aircraft  voltage  will  not  effect  instrument  calibrations  with  this 
arrangement.    However,  it  was  desired  to  improve  the  sensitivity  of  the 
angle  of  attack  and  normal  acceleration  measurements;  therefore,  28 
volts  were  employed  as  full  scale.     Since  the  telemeter  will  accept  only 
zero  to  five  volt  signals,  a  bias  battery  was  employed  to  set  the  steady 
state  potentiometer  inputs  within  this  range.     This  has  the  disadvantage 
that  a  change  in  the  28  volt  bus  will  change  the  steady  state  output  of  the 
angle  of  attack  vane  and  the  normal  accelerometer*    However,  the  sensi- 
tivity will  not  be  effected. 
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Ground  Station  Equipment 

The  ground  station  equipment  included  a  telemeter  receiver,  a 
magnetic  tape  recorder,  an  analog  computer  and  a  Sanborn  model  record- 
er.    These  components  are  shown  in  Figs.  4  and  5. 

The  telemeter  receiver  was  built  by  the  Applied  Science  Corpor- 
ation of  Princeton.     This  unit  is  designed  to  receive  the  coded  output  of 
the  telemeter  transmitter,  decode  the  signals,  mix  channels  ae  desired 
and  transmit  the  data  to  a  recording  device. 

An  Ampex  Model  309C  dual  track  tape  recorder  was  employed  to 
record  the  data.    This  provided  a  convenient  means  of  retaining  data  until 
ready  for  further  analysis.     The  fact  that  the  data  was  in  an  electrical 
form  was  a  great  advantage  when  the  use  of  an  analog  computer  is  contem- 
plated.   All  voice  transmissions  were  also  recorded  on  the  second  tape 
channel. 

A  Goodyear  Aircraft  Corporation  Model  L  3   (GEDA)  linear 
electronic  differential  analyser  was  employed  to  analyse  the  data.     This 
analog  computer  provides  twenty -four  stabilised  DC  amplifiers  with  an 
open  loop  gain  of  5  x  10.      The  computer  had  a  guaranteed  accuracy  of 
one  percent.     Provisions  for  accurately  setting  the  computer  board  poten- 
tiometers was  available. 

A  Sanborn  Model    154-100  B  four  channel  recorder  was  used 


10 


for  visual  presentation  of  the  data.     This  instrument  had  a  high  natural 
frequency,   and  the  low  frequency  response  was  flat  to  zero. 

Instrument   Calibration 

Three  components  of  the  instrumentation  system  required  cali- 
bration.   These  were  the  normal  accelerometer,   pitch  rate  gyro  and  angle 
of  attack  vane.     The  aircraft  airspeed  indication  and  altimeter  were  not 
calibrated  as  previous  investigation  indicated  small  errors  in  these 
systems. 

The  accelerometer  was  calibrated  by  tilting  it  through     ■*  90°. 
The  angle  of  tilt  and  the   percent  of  full  scale  voltage  were  recorded. 
This  checked  very  closely  with  the  factory  calibration  which  was  there- 
fore accepted  as  valid.     Fig.  7  shows  the  sine  of  the  angle  of  tilt  versus 
the  potentiometer  output  in  percent  full  scale  voltage.     The  slope  of  this 
curve  was  found  to  be  2.  115   g's  for  full  scale  voltage.     Since  the  instru- 
ment was  canted  at  4b°,   this  value  was  divided  by  the  sine  of  45°  to 
give  2.995  g's  for  full  scale.     As  installed  in  the  aircraft  28  volts  were 
applied  across  the  accelerometer  potentiometer  with  a  five  volt  portion 
of  this  representing  the  telemetering  full  scale.     As  received  at  the  tele- 
metering ground  station,   full  scale  was  therefore:      5/28  x  2.995    * 
.534  g»s. 

Part  way  through  the  flight  test  program,   it  was  discovered 
that  when  releasing  the  bomb  electrically,   the  bomb  release  solenoid 
reduced  the  battery  bus  voltage  by  .6  volts*     This  reduced  voltage  lasted 
as  long  as  the  spring  loaded  bomb  release  switch  was  held  down.      This 
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effectively  reduced  the  accelerometer  zero  by  12%  of  telemetering  full 
scale.  This  is  indicated  by  the  negative  peak  on  the  flight  4  accelerom- 
eter trace  in  Fig.    14. 

Henceforth,   the  bomb  was  released  manually.     However,   in  the 
two  additional  drops  which  were  accomplished  the  accelerometer  was 
found  to  have  a  dead  zone  in  the  steady  state  position.     Time  would  not 
permit  additional  drops.     Careful  recalibration  of  the  instrument  showed 
it  was  linear  outside  the  dead  zone,   and  the  zero  was  displaced  .  6  g. 
The  new  calibration  is  shown  on  Fig.    8  as  .518  g's  for  telemeter  full 
scale. 

The  dynamic  characteristics  of  the  accelerometer  were  deter- 
mined as  follows.     The  instrument  was  placed  on  a  board  with  the  sensi- 
tive axis  vertical.     The  support  was  then  sharply  removed,   and  the  tran- 
sient response  of  the  instrument  recorded  during  its  free  fall.     A  drop  of 
approximately  four  inches  was  sufficient  for  the  instrument  to  reach  a 
zero  g   steady  state.     The  instrument  dynamic  characteristics  were  de- 
termined by  matching  an  analog  computer  trace  to  the  instrument   re- 
sponse.    The  damping  ratio  of  .481  and  the  natural  frequency  of  6.00 
cycles  per  second  could  then  be  found  from  the  computer  potentiometer 
settings. 

The  pitch  rate  gyro  was  calibrated  by  placing  it  on  a  constant 
speed  turntable  and  by  recording  the  revolutions  per  minute  and  the  poten- 
tiometer output.     In  order  to  decrease  the  full  scale  range,   a  portion  at 
each  end  of  the  potentiometer  was  shorted  with  metallic  paint.     The  final 
calibration  curve  is  shown  on  Fig.   9.     Full  scale  deflection  represents 
*•   8.  5°  per  second. 
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The  gyro  dynamic  characteristics  were  determined  by  deflect- 
ing the  rotor  and  by  recording  its  return  to  the  equilibrium  condition. 
The  analog  computer  was  used  to  determine  its  damping  ratio  of  «0597 
and  its  natural  frequency  of  4.  58  cycles  per  second. 

The  angle  of  attack  vane  was  calibrated  in  two  ways.      On  the 
ground,   a  protractor  was  mounted  in  such  a  way  that  the  angular  deflec- 
tion of  the  vane  could  be  read.     The  angle  and  potentiometer  output  were 
recorded  and  are  shown  on  Fig.    10.     An  in-flight  calibration  was  made 
by  uoing  the  telemetering  receiver  to  record  the  vane  position.     With 
the  aircraft  in  steady  level  flight,   the  aircraft  attitude  was  read  with  a 
bubble  level.     By  flying  at  varying  speeds,   a  direct  plot  of  angle  of 
attack  versus  percent  telemeter  full  scale  is  then  available.     Fig.    10 
shows  both  calibrations.     The  ground  calibration  indicated  telemetering 
full  scale  was  9.  04  degrees  while  the  inflight  calibration  gave  7.95  de- 
grees for  telemetering  full  scale.     The  dynamic  characteristics  of  the 
instrument  were  not  determined  due  to  poor  data  obtained  from   this 
instrument. 

In  addition  to  the  instrument  calibration,   it  was  necessary  to 
determine  the  aircraft  weight,   e.g.  location,  and  moment  of  inertia. 
Wheel  scales  were  employed  to  determine  the  aircraft  total   weight. 
The  aircraft  was  then  placed  on  two  wing  jacks  and  supported  at  the  tail 
with  an  accurate  balance  scale.     The  e.g.   location  was  then  determined 
from  the  location  of  the  support  points.     These  measurements  were 
made  with  the  pilot  in  the  aircraft  and  the  wheels  raised.     The  test  air- 
craft was  equipped  with  a  55  pound  sliding  lead  weight  which  was 
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mounted  in  a  hollow  tube  in  the  fuselage  and  which  could  be  moved   17.2 
feet  through  a  crank  arrangement  in  the  cockpit.     The  flight  tests  were 
made  with  the  lead  weight  in  the  aft  and  in  a  mid  position.     The  center  of 
gravity  under  these  conditions  is  29.2  and  32.5%   mac.     While  measuring 
horizontal  distances  with  the  aircraft  in  the  flight  attitude,   the  bomb  c«  g. 
by  3.4  and  1. 1  inches  for  the  aft  and  mid  e.g.   positions  respectively, 
Because  the  bomb  c.  g.  is  at  a  reasonable  distance  vertically  from  the 
aircraft  e.g.,   these  distances  are  dependent  upon  an  accurate  attitude 
measurement  of  the  flight  conditions  at  the  time  of  the  bomb  drop.     This 
was  accomplished  by  a  flight  with  the  bomb  installed  using  a  bubble  level. 
The  aircraft  was  then  placed  in  the  same  attitude  while  on  jacks  on  the 
ground  and  the  aircraft  and  bomb  e.g.   positions  were  determined.      The 
moment  of  inertia  of  the  aircraft  v/as  determined  by  oscillating  the  air- 
craft about  the  jack  points  with  a  spring  mounted  on  the  tail.      Fig.    11 
shows  the  configuration  employed.     The  damping  under  these  conditions 
proved  to  be  negligible.     The  natural  frequency  was  determined  by  timing 
50  cycles  with  a  stop  watch.     Knowing  the  natural  frequency  and  the 
spring  constant  of  the  system,   the  moment  of  inertia  about  the  jack  points 
was  determined.     Calibration  of  the  springs  showed  them  to  be  linear  for 
the  load  range  required.     The  moment  of  inertia  was  then  transferred  to 
the  aircraft  Co  g.      The  moments  of  inertia  about  the  aircraft  e.g.  were 
found  to  be  3536  and  3126  slug  ft.2  with  the  e.g.   in  the  aft  and   mid 
positions  respectively. 
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PROCEDURE 

The  procedures  employed  in  this  investigation  were  basically- 
very  simple.     A  steady  state  flight  condition  was  established;  the  con- 
trols were  locked;   and  the  bomb  was  dropped.     The  initial  airspeed, 
pressure  altitude,   and  temperature  were  recorded  by  the  pilot.      The 
transient  response  of  the  aircraft,   as  sensed  by  the  instruments,   was 
transmitted  through  the  telemetering  units  and  recorded  on  magnetic 
tape. 

In  detail,   the  procedures  were  carried  out  as  follows.     Since 
it  was  desired  to  recover  the  bomb,   a  water  drop  was  impractical.     In 
order  to  drop  over  land  with  reasonable  safety,   only  low  altitude  drops 
at  less  than  500  feet  were  considered  practical.     In  addition,   a  high  alti- 
tude drop  would  require  a  major  excavation  to  recover  the  bomb. 

Since  very  smooth  air  was  known  to  be  absolutely  necessary, 
it  was  decided  to  make  all  drop  flights  very  near  sunrise.     With  gener- 
ally good  weather,  the  air  was  smooth  below  500  feet  at  approximately 
every  other  day.     This  gave  an  average  of  about  two  satisfactory  morn- 
ings a  week  to  attempt  the  drop.     Winds  of  three  to  five  miles  per  hour 
rendered  conditions  unsatisfactory  for  the  low  altitude  test  flight. 

The  initial  intent  was  to  secure  good  bomb  drop  data  for  three 
c.  g.   positions  at  a  single  airspeed.     In  addition,   it  was  intended  to 
record  the  aircraft  response  to  a  small  elevator  step  for  at  least  one 
set  of  conditions  common  with  a  bomb  drop.     A  total  of  nine  bomb  drops 
were  made.     Due  to  instrumentation  and  weather  problems,   no  com- 
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pletely  satisfactory  drop  was  accomplished.     However,   two  drops  were 
made  where  it  was  felt  that  the  unknowns  could  be  removed  with  reason- 
able confidence.     Also,   one  elevator  step  was  recorded  that  was  felt  to 
be  usable.     The  data  for  these  runs  are  shown  in  Figs.    14  through   16. 

The  pre-flight  checks  required  for  a  test  flight  were  practically 
non-existent.     It  was  found  desirable  to  check  that  the  aircraft  and  instru- 
mentation were  in  good  condition  the  afternoon  prior  to  a  test  flight.    The 
weather  predictions  were  checked  and  arrangements  made  to  gain  access 
to  all  necessary  equipment  for  the  drop  on  the  following  morning. 

Prior  to  take-off,    radio  contact  was  made  with  the  ground  sta- 
tion to  ascertain  proper  functioning  of  the  telemetering  installation. 
After  take-off,   the  lead  weight  was  placed  in  the  desired  location,  and 
the  aircraft  was  trimmed  for  level  flight.     When  proper  power  settings 
were  established  the  ground  station  was  advised  that  a  run  was  being 
commenced.     The  elevator  was  locked  approximately  30  seconds  prior 
to  the  bomb  release  with  a  clamp  placed  on  the  control  column.     If  air- 
speed and  attitude  remained  constant  for  the  final  30  seconds  and  no 
gusts  were  encountered,   the  bomb  was  released.     It  was  found  that  200 
to  400  feet  was  the  most  practical  altitude  range  for  bomb  release. 
Higher  altitudes  made  bombing  accuracy  uncertain.     Lower  altitudes 
were  found  to  lack  the  necessary  smooth  air  and  required  excessive 
care  in  maintaining  safe  clearances. 

Once  take-off  was  made,   the  bomb  was  dropped  regardless  of 
any  instrumentation  or  telemetering  problems.     This  required  that 
streamlining  be  remounted  on  the  bomb,   but  it  was  felt  that  it  was 
unwise  to  make  a  landing  with  the  bomb  installed. 
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DISCUSSION  AND  RESULTS 

Theory 

In  studying  Equations  of  Motion  of  the  longitudinal  transient 
response  of  an  aircraft,   it  is  conventional  to  assume  that  velocity 
remains  constant.     Also,  the  aerodynamic  forces  can  be  expressed  by 
a  Taylor  expansion  about  the  equilibrium  condition  retaining  only  linear 
terms.     Under  these  conditions  Ref.  1    develops  the  equations  of  motion 
as  follows: 

(1)  ta       a    +    da     •  dG   *    i l  (t) 

(Z)       Cma    a    ♦    Cm^    da    ♦   Cmde   dO    -    hd2G   «   f2(t) 

In  addition,   normal  acceleration  can  be  expressed  as: 


*3*      _-2.   n   *     dG   -  da 


The  particular  forcing  function  in  this  investigation  can  lead 
to  some  confusion  in  considering  initial  conditions  and  in  placing  the 
equations  in  non-dimensional  form.     In  this  investigation,   the  aircraft 
transient  response  was  induced  by  dropping  a  weight  with  the  controls 
fixed.     Clearly,   before  the  weight  is  dropped,   the  lift  is  equal   to  the 
total  weight  of  the  aircraft  and  bomb.     However,   the  mass  of  the  air- 
craft during  the  transient  response  does  not  include  the  bomb  weight. 
To  avoid  confusion,   it  is  best  to  think  of  the  aircraft  configuration  for 
t  <  0    to  be  that  of  the  basic  aircraft  (no  bomb)  with  its  c.  g„   deter- 
mined,   plus  the  effects  of  the  bomb  acting  in  relation  to  the  basic  air- 
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craft  c.  g.     I€  the  airplane  plus  bomb  configuration  iB  trimmed  at    t  <  0, 
then  releasing  the  bomb  ie  equivalent  to  applying  an  up  force  at  the  bomb 
c,  go      Thus,   one  forcing  function  in  the  lift  equation  is  an  up  force  equal 
to  the  bomb  weight,   and  in  the  moment  equation,   we  have  a  positive 
pitching  moment  equal  to  the  bomb  weight  times  the  distance  between 
the  basic  aircraft  c.  go  and  the  bomb  c„  g0 

b    s    bomb  weight 


(4) 

o    »                    b 

Wb 

(5) 

psvz 

2  Wb 
pSV* 

(6) 

Let:         CL    1      » 
o 

2  W 

pSV2 

It  is  pointed  out  that  the  initial  angle  of  attack  is  prescribed 
by     Cj^       as  defined  above.     However,   the  aerodynamic  time    (t)  and 
the  dimensionles8  moment  of  inertia  (h)  are  defined  by  using  the  air- 
craft mass  during  the  transient  response.     It  is  impossible  therefore 
to  duplicate  all  conditions  using  any  forcing  function  other  than  that 
described  previously. 

In  addition,   it  is  necessary  to  examine  the  expression  for 
normal  acceleration.     This  is  derived  from  the  expressions  developed 
in  Ref ,    1  : 


(7)        F      c    nW    «    mV(G-fl) 
z 


If  this  expression  is  changed  to  non-dimensional  form,   the 
following  results: 
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(8)       tdW    j    mV(dG-da) 
T 


(9)       n      w   ■     *    (de  -  da) 


pSV 
(10)     n  C^o       -     (dG  -  do) 


The  forcing  functions  in  dimensional  form  are  the  bomb 
weight  in  the  lift  equation;  the  bomb  weight  times  the  distance  between 
the  aircraft  and  the  bomb  center  of  gravity  in  the  moment  equation. 
In  non-dimensional  form,  'the  forcing  functions  become: 

For    t  <  to  fj(t)  *   f^t  ■  0 

For    t  >  to 

w  w        C,    1 

(ID     MO    «    -_b         *    -    _%         l-0 

Tsv1  w     "~"r~ 


o 


(12)     f2(t)    -     2  Wb*b      ■    -Wb       CL   1      i^ 


Define  the  quantities: 


Kn  W,        CT     1  Kn        m      W,        C.     1      L 

AC.     »        b  L        ,      AC        «         b  L,  d 

L         -^  o  m  -^  o      — 


Then:        f.  (t)    *    -ACl,        and      f,  (t)    «    -  ACm 
i  a  c,  m 


All  time  histories  obtained  in  transient  responses  are  in 
real  time 


19 


(13) 

a 

2  T 

a 

'            *                 AC 

+  a   -   e  «    -        L 

T 

(14) 

c 

2 

T 

a 

c                   c 

T                                  T 

• 

e 

•  • 

-  h  0    = 

AC 
m 

T2 

These  expressions  are  identical  to  those  obtained  using  an  elevator 
step  input  except  for  the  definitions  of  the  forcing  functions.     For  an 
elevator  step  input: 

AC  C  .  AC.  AC  n 

m    a        m.       0      ,  JL    »  m        C 

It  is  pointed  out  again  that  all  conditions  cannot  be  matched  for  the 
two  types  of  forcing  functions.     In  practical  term,    initial  angle  of 
attack  and  velocity  cannot  be  identical  for  the  two  cases  at  the  same 
density  altitude „ 

There  is  an  additional  point  of  view  which  can  be  taken  about 

the  initial  conditions  of  the  bomb  drop.     Equations  (1)  and  (2)  are 

perturbation  equation  So     By  considering    f .  (t)    and    f.>  (t)    zero  before 

the  bomb  drop,   this  establishes    a      and    dG      as  zero.      In  a  linear 
r  o  o 

perturbation  problem  where  the  initial  and  final  solutions  are  both 
steady  state  conditions,   it  is  arbitrary  as  to  which  is  considered  the 
undisturbed  state.     It  is  equally  valid,   in  the  aircraft  transient  re- 
sponse to  a  bomb  drop,   to  consider  the  final  steady  state  as  the  undis- 
turbed condition  with    a    and     dG     both  zero.       The  physical  motion 
of  the  aircraft  is  a  steady  pitch  rate. 
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This  can  be  more  clearly  demonstrated  with  the  use  of  equa 
tions.     The  variables  can  be  expressed  as  the  sum  of  a  steady  state 
and  a  transient  solution. 


a    »   a        +   a.  dGsde+dG. 

86  t  SS  t 


Under  the  previous  arrangement  when    t    =    O: 


a      «    O      a       *   -a.  dG     *    O       dG       =    -dG, 

o  ss  t  o  ss  t 


If  the  final  steady  state  is-- defined  as  the  undisturbed  condition,   the 
equations  become  : 


ase    *    dess    *    °  a    *    at  dc?    s    det 


Under  these  conditions    f.  (t)    and    i.,  (t)     are  zero  after  the  bomb  drop 
and  are  the  negatives  of  the  previous  forcing  functions  before  the  bomb 
drop.     That  is,   the  aircraft  before  the  bomb  drop  has  a  force  and 
moment  acting  until  release.     The  initial  perturbation  angle  of  attack 
and  pitch  rate  are  unknown  in  this  case.     An  analytic  solution  of  the 
equations  will  show  that  the  steady  state  solutions  in  the  first  case 
are  the  negatives  of  the  initial  condition  in  the  second  case. 

The  equations  for  the  steady  state  solutions  can  be  found  from 
equations   (1)  and  (2)  for  the  original  arrangement. 

(15)    de.      =     "La/z     V*>    -Cma    fi(t) 

ss      — c — :  c,   / — c 

ma    *      Vz       mdb 


(16)    aBS    -     V"     ♦     °mde    fl(t) 
"C IT, — /     C 


ma    +     \L       md£ 
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The  initial  values  for     d  9     and       a      in  the  second  case  are  repre- 
sented by  identical  expressions  as  above.    However,   the  forcing 
functions  change  signs. 

In  the  real  aircraft  problem,  there  is  some  disadvantage  in 
considering  the  final  state  as  the  undisturbed  condition.     The  aircraft 
may  not  reach  a  complete  steady  state  in  the  time  period  where  veloc- 
ity is  essentially  unchanged.     There  is  also  a  mental  difficulty  in 
thinking  of  a  steady  pitch  rate  as  an  undisturbed  condition  in  relation 
to  level  flight. 

.Linear   Dependence 

The  determination  of  unique  values  for  stability  derivatives 
from  flight  test  data  is  a  problem  that  arises  in  any  experimental  pro- 
gram.    The  stability  derivatives  which  can  be  determined  uniquely  is 
controlled  by  the  linear  dependence  of  the  equations  describing  the 
system «     This  subject  is  considered  in  detail  in  Refs.  Z  and  3.     For 
the  longitudinal  motions  considered  in  this  report,   the  problem  can 
be  resolved  without  great  difficulty. 

In  equations  (13)  and  (14),   it  is  considered  that  all  flight 

variables  and  the  forcing  functions  are  known  and  are  not  zero. 

Under  these  circumstances,   there  is  no  linear  dependency  between 

the  two  equations.    All  stability  derivatives  and  the  dimensionless 

moment  of  inertia  (h)   can  theoretically  be  determinedo     The  ability 

to  determine     C  and     C  independently  by  dropping  a 

m  ._  m  ,  ^ 

da  d6 

weight  is  the  basic   reason  this  investigation  was  undertaken.      The 
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use  of  the  elevator  as  a  forcing  function  does  not  accomplish  this  for 
one  reason.     The  tail  force  introduced  into  the  lift  equation  is  so  small 
that  its  effect  on  the  aircraft  motion  is  negligible.     Equation  (13)  may 
effectively  be  equated  to  zero  and  a  relation  between  the  variables   a.  , 
a    and     0      results.     Therefore,  the  variables  in  equation  (14)  are  not 
linearly  independent,  and  it  is  impossible  to  find  all  the  stability  deriv- 
atives in  this  equation  uniquely.    A  complete  development  of  what  can 
be  found  is  contained  in  Ref.  2. 

In  using  a  bomb  drop  as  a  forcing  function,  there  is  an  addi- 
tional advantage.     The  exact  nature  of  the  forcing  functions  are  known 
as  a  force  and  a  moment.     With  an  elevator  step,  this  is  not  the  case. 
It  is  pointed  out  that  the  moment  produced  by  the  bomb  release  does  not 
have  to  be  known  to  solve  the  problem  if   h    is  known. 

The  facts  stated  above  can  be  demonstrated  by  the  use  of  the 
equations  a6  follows.     If  equation  (13)  is  used  to  eliminate     0      from 
equation  (14),  the  following  results: 


(c-« 


(17)    \     "a    T  ld0j  o?    +       mda  +    md0  -h    La 


♦ 

LQJ             md0 

2 

/ACl    C™d0 

T 

/           , 

» 

a 


)-<»(r^]-^? 
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If    A  CL     and    ^  ^m     are  known  and  are  not  zero,  all  quantities  in 
parenthesis  can  be  determined  with  a  time  history  of    a  ,    a    and  a.  . 
If     h     is  known,   it  is  not  necessary  to  know  the  value  of     ^^m     to 
determine  these  quantities.     Inspection  of  the  quantities  in  parenthesis 
shows  that  all  derivatives  are  determined  uniquely.     However,  if 
^CL    »    0,  as  is  the  case  for  an  elevator  input,  the  derivatives  cannot 

A   C 

be  found  uniquely  even  if     h     and  m     are  both  known.     It  is  also 

pointed  out  that     ^cm     can  be  identically  zero  in  the  case  of  the  bomb 
drop  without  destroying  the  uniqueness  of  the  solution. 

Preliminary  Investigation 

The  previous  analysis  verified  that  a  solution  was  possible  for 
all  longitudinal  stability  derivatives  by  dropping  a  weight.    However,  it 
was  still  necessary  to  verify  that  the  transient  response  of  the  aircraft 
would  be  of  sufficient  magnitude  that  the  instrumentation  available  could 
sense  the  flight  variables  with  the  necessary  accuracy.     It  was  deter- 
mined that  about   500  pounds  could  be  carried  externally  on  the  te6t  air- 
craft.    This  represented  approximately  20%  of  the  basic  aircraft  weight. 
It  was  estimated  that  the  bomb  and  aircraft  centers  of  gravity  would  be 
approximately  eight  inches  apart.     Refs.  4  and  5  provided  reasonable 
estimates  for  the  longitudinal  stability  derivatives.     "With  this  informa- 
tion,  it  was  possible  to  estimate  the  response  of  the  aircraft  when  the 
bomb  was  dropped. 
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Analog  of  the  Test  Aircraft 

It  was  decided  that  the  use  of  an  analog  computer  was  the  most 
convenient  method  to  make  a  preliminary  analysis.     Equations  (13)  and 
(14)   were  set  up  on  the  computer  and  the  time  history  of  the  aircraft 
simulated  flight  parameters  were  placed  on  a  Sanborn  recorder.     In 
addition,  equation  (7)  was  set  up  to  give  normal  acceleration.    See  Fig. 
12  for  a  schematic  wiring  diagram  of  the  analog  computer  setup.     This 
arrangement  allowed  variation  of  the  forcing  functions  and  all  stability 
derivatives  as  desired.     This  arrangement  assumed  the  aircraft  was 
perfectly  described  by  the  equations  and  that  the  instruments  would  per- 
fectly  record  the  flight  variables    a  ,    9     and    n. 

Solution  for  Stability  Derivatives  from  Analog  Data 

It  was  found  that  the  stability  derivatives  could  be  found  with 
reasonably  good  accuracy  using  the  very  simple  statistical  methods  dis- 
cussed in  Ref.  6.    However,  it  was  found  that  greater  accuracy  in  much 
less  time  could  be  achieved. by  curve  matching  two  computer  solutions. 
One  person  would  set  in  arbitrary  stability  derivatives  and  produce  a 
time  history  of  aircraft  motions.     A  second  person  would  then  reproduce 
the  curves  as  closely  as  possible  and  determine  the  stability  derivatives 

from  the  potentiometer  settings .    Excellent  results  were  obtained  in  this 

C  C 

manner  for  the  derivative       mdQ  .      Good  results  were  found  for        ma 

c  c 

while  only  fair  accuracy  could  be  expected  for        nida     and        Lq;  • 
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It  was  al6o  discovered  that  the  variation  in  angle  of  attack 

would  be  quite  small.     This  is  of  course  related  to  the  difficulties  in 

Cm  Ct 

determining  a.      and  a  .      Improved  results  could  be  obtained 


by  increasing  m,     that  is,  moving  the  bomb  further  from  the  air- 

craft center  of  gravity.  Unfortunately,  the  location  of  the  bomb  on  the 
aircraft  was  quite  closely  fixed  by  the  landing  gear  wheel  wells.  The 
remaining  solution  was  to  shift  the  aircraft  center  of  gravity.  This  is 
not  nearly  as  satisfactory  since,  in  practice,   shifting  weight  aft  in  the 

Q 

aircraft  increases  its  moment  of  inertia  and  rapidly  decreases       i*^  . 
Both  of  these  effects  tend  to  reduce  the  size  of  the  aircraft's  response. 
The  analog  computer  demonstrated  that  moving  the  aircraft  e.g.  aft 
increased  the  response  qualitatively  about  half  as  much  as  moving  the 
bomb  c,  g«  forward.     The  character  of  the  changes  were  distinct  how- 
ever. 


In  addition,  the  use  of  elevator  input  to  create  a  sizeable 
moment  input  was  investigated.    It  was  found  in  this  case  that  the  re- 
sponse  was  effected  most  strongly  by        Lq,      and        mfl  .        It  was 
known  from  the  previous  theoretical  discussion  that  no  unique  solution 
existed  for     C™a ,     Cmda     and     CmdG  •       However ,  the  computer 
demonstrated  that  if  a  reasonable  approximation  was  available  for  rate 

Q 

derivatives,  good  results  could  be  obtained  for        mG!      as  well  as 
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c  c  c 

Using  both  type  inputs,       mae  ,        rafi     and        La      could  be 

c 
reproduced  to  a  close  percentage  of  the  initial  value.         Tnda      could 

not  be  reproduced  with  a  high  percentage  of  accuracy,  but  in  terms  of 

absolute  value,  it  was  well  defined.     This  difficulty  with        mdQ!      is 

somewhat  inherent  due  to  its  small  size  compared  to        m^a        and 

^Tna  .      As  Ref.  7  states,  the  accuracy  with  which  a  stability  derivative 

can  be  found  is  determined  by  the  magnitude  of  its  influence  on  the  motion 

Q 

of  the  aircraft.  Tndd      k&s  a  relatively  small  influence  on  the  air- 

craft's motions  and  therefore  will  not  be  found  with  high  accuracy. 


The  variables  which  were  to  be  measured  included  normal 
acceleration,  pitch  rate  and  angle  of  attack.    Any  two  of  the  three  will 
provide  a  complete  solution  a6  is  shown  by  equations  (1),  (2),  and  (3). 
Under  the  conditions  6et  into  the  analog  computer  for  the  bomb  drop, 
the  time  histories  indicated  the  following  change  from  the  steady  state 
in  the  variables  to  be  measured: 

Normal  acceleration:      0  to   .  4  g 

Pitch  rate:  0  to  6°  per  second 

Angle  of  attack:  -.  5°  to    +2° 

It  was  felt  that  acceleration  and  pitch  rate  could  be  measured 
in  this  range  with  reasonable  accuracy.      However,  it  was  felt  that 
angle  of  attack  could  probably  not  be  measured  with  sufficient  accuracy 
to  give  any  useful  information. 
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Results  of  Preliminary  Investigation 

The  following  facts  were  indicated  by  the  use  of  the  analog  com- 
puter in  the  preliminary  investigation. 

a)  The  longitudinal  stability  derivatives  could  be  found  uniquely  by 
dropping  a  weight  approximately  20%  of  the  aircraft  weight. 

b)  The  weight  should  be  located  further  from  the  aircraft  center  of 

gravity.    A  distance  of   1.00  to  1.75  feet  would  be  satisfactory  for  the 

Navion.       (i    «    .176  to.  305) 
c 

Q 

c)  The  derivative        mdQ       would  be  determined  most  accurately. 
mda      would  be  found  with  only  fair  percentage  accuracy  but  good  abso- 
lute accuracy.     The  use  of  an  elevator  step  would  assist  in  determining 

Tna      ttd        J->a  '       However,  as  mentioned  previously,  all  conditions 
of  flight  cannot  be  duplicated  for  the  two  forcing  functions. 

d)  It  was  likely  that  no  reliable  angle  of  attack  information  would 
be  obtained. 

e)  The  best  means  of  obtaining  a  solution  was  probably  matching 
analog  data  with  the  measured  response  of  the  aircraft. 

The  general  technique  for  finding  a  set  of  stability  derivatives 
by  curve  matching  was  believed  to  be  as  follows: 

a)      With  the  best  available  estimates,  create  time  histories  with 
the  analog  computer  of  all  flight  variables  which  are  available  from  the 
test  flight.     The  step  forcing  functions  are  easily  created  by  using  the 
reference  voltage  of  the  computer. 
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b)  Vary  each  stability  derivative  separately  and  make  an  overlay 
for  each  derivative  showing  its  effect  on  the  time  histories. 

c)  Carry  out  this  procedure  for  each  type  of  forcing  function 
considered. 

d)  Use  the  bomb  drop  time  history  and  create  the  best    match 
possible. 

e)  Shift  to  an  elevator  step.     Leave  the  derivatives        md0      and 
"Ma      fixed,  and  create  the  best  match  possible. 

f)  Repeat  steps  (d)  and  (e)   until  the  best  possible  match  is 
achieved. 

Some  small  allowance  could  be  made  for  the  difference  in  flight 
conditions  between  the  two  forcing  functions.    If  identical  indicated  air 
speeds  are  used,  the  initial  angle  of  attack  for  the  elevator  step  will  be 
approximately  two  degrees  less  than  that  for  the  bomb  drop.     This  con- 
dition is  largely  relieved  by  the  fact  that  the  angle  of  attack  increases 
three  to  four  degrees  during  the  elevator  step  while  it  shows  very  small 
change  during  the  bomb  drop. 

Results  Using  Flight  Test  Data 

As  previously  mentioned  in  the  report,  no  completely  satisfac  - 
tory  data  was  obtained.     Two  bomb  drops  and  one  elevator  step  were 
selected  as  providing  usable  information.    Figs.   14  through   16  show 
this  data.     The  flight  conditions  for  these  runs  are  recorded  in  Table  U. 
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When  the  bomb  rack  installation  was  completed,  it  was  found 
that  the  aircraft  and  bomb  centers  of  gravity  were  considerably  closer 
together  than  had  been  initially  estimated.    This  had  a  twofold  disadvan- 
tage.    As  previously  discussed,  it  reduces  the  size  of  the  aircraft 
response,  making  the  instrumentation  requirements  more  severe. 
Also  it  makes  the  determination  of  the  size  of  the  moment  forcing  func- 
tion more  difficult.     The  horizontal  distance  between  the  aircraft  and 
bomb  centers  of  gravity  depend  quite  strongly  on  the  aircraft  attitude 
when  the  centers  of  gravity  are  very  close  together  horizontally  in  level 
flight.    A  one  degree  error  in  aircraft  attitude  created  approximately  a 
20%  error  in  the  size  of  the  pitching  moment  when  the  aircraft  e.g. 
was  completely  aft.     The  situation  was  even  worse  as  the  aircraft  cen- 
ter of  gravity  moved  forward.    Since  no  instrument  recorded  flight  atti- 
tude, it  was  felt  that  no  accurate  moment  forcing  function  could  be 
determined.    Measurement  of  aircraft  attitude  with  the  bomb  installed 
showed  a  variation  too  large  to  be  ignored.    It  was  therefore  decided  to 
regard  the  moment  input  as  an  unknown.    As  brought  out  in  the  theory, 
the  problem  6till  has  a  unique  solution;  however  practical  application 
showed  it  was  much  more  difficult. 

The  streamlining  on  the  bomb  was  initially  installed  to  reduce 
the  aerodynamic  drag.     This  was  thought  to  be  desirable  in  order  to 
estimate  the  pitching  moment  more  accurately.    An  estimate  of  .  5  for 
the  drag  coefficient  of  the  streamlined  bomb  was  made  based  on  its 
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cross  sectional  area.     This  estimate  was  taken  from  Ref.   1.      The 
moment  resulting  from  the  bomb  drag  was  determined  to  be  quite  small 
and  therefore  unimportant  since  the  pitching  moment  was  not  accurately 
known  anyway.     The  streamlining  was  retained  throughout  this  investi- 
gation in  an  effort  to  reproduce  results.    It  is  not  considered  to  be 
necessary  for  any  future  investigation  of  this  nature  using  low  speed  air 
craft.    It  should  be  pointed  out  that  it  has  been  assumed  that  there  is  no 
aerodynamic  effect  on  the  aircraft  due  to  the  bomb's  presence. 

The  matching  of  the  time  histories  for  two  bomb  drops  and  one 
elevator  step  are  shown  in  Figs.  17  through  19.  The  stability  deriva- 
tives obtained  through  these  curves  are  recorded  in  Table  III. 

In  the  analog  computer  solution,  the  dynamic  characteristics 
of  the  accelerometer  and  pitch  rate  gyro  were  included.     The  computer 
circuit  diagrams  for  these  instruments  are  shown  in  Fig.    13.     The  use 
of  filters  to  remove  noise  from  the  test  data  wa6  considered.     In  gen- 
eral, it  was  found  that  the  stray  noise  could  be  removed  by  eye  better 
than  by  filtering  when  curve  matching  was  employed  as  a  method  of 
solution. 

The  effect  of  each  stability  derivative  on  the  computer  time 
histories  is  6hown  in  Figs.  20  through  24.     It  can  be  seen  that    ^^m, 
^La      and        m#      are  not  clearly  defined  by  the  use  of  the  bomb  drop. 
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However,  d9       and  da      appear  to  be  reasonably  well  defined. 

C  C 

The  elevator  step  improve©  the        Lce      and        Tna  ,    but  the  failure  to 

reproduce  flight  conditions  prevents  high  accuracy.     Figs.  25  through 
29  show  the  effect  of  each  stability  derivative  on  the  response  to  an 
elevator  step. 

In  spite  of  the  instrumentation  difficulties  previously  discussed, 
the  numerical  results  of  this  investigation  are  considered  to  be  reason- 
ably good.     In  the  elevator  step  shown  in  Fig.    19,  there  are  no  ambigui- 
ties in  angle  of  attack  or  pitch  rate.     Table  IV  records  the  derivatives 
as  established  in  this  investigation  and  those  arrived  at  in  Refs.  4  and  5. 

c  c  c 

The  derivatives        md0  ,         mda     and        ^a      established  in  this  inves- 
tigation fall  within  the  band  of  previous  evaluations.     Conditions  were 
not  reproduced  sufficiently  to  allow  a  good  comparison  for  mQ»  . 

The  indications  are  that        rna     as  established  in  this  report  is  at 

C— 
least   10%  high.  Tna     was  determined  with  the  aircraft  e.g.  located 

at  29-2%  and  32.5%  mac.     This  gives  a  theoretical  decrease  of  .  17  in 

the  magnitude  of        m(¥  .       The  experimental  decrease  was  .22. 

The  calculated  and  experimental  values  of  ^  cm  for  the 
bomb  drops  agreed  within  15%.  In  both  cases,  the  calculated  values 
were  greater.      See  Table  III. 
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Proposed  Improvements 

The  basic  improvements  which  should  be  made  have  already- 
been  discussed.    It  is  an  obvious  fact  that  the  instruments  employed 
should  be  reliable  and  cover  the  desired  range  without  modification.    If 
telemetering  is  employed,  it  is  a  great  advantage  to  use  the  telemetering 
full  scale  voltage  as  the  reference  voltage  for  all  potentiometer  trans  - 
due ere. 

It  would  be  of  considerable  assistance  if  aircraft  attitude  and 
elevator  position  were  also  measured.     While  neither  is  required,  they 
both  help  to  determine  the  quality  of  the  steady  state  established  prior  to 
the  bomb  release.    In  addition,  aircraft  attitude  could  be  used  to  estab- 
lish the  moment  forcing  function.    Elevator  position  would  verify  that  the 
controls  remained  fixed  during  the  bomb  release  and  could  be  used  as  a 
forcing  function  if  an  elevator  input  was  employed.     This  would  also 
allow  the  determination  of        m  6      rather  than     ^  ^m     for  the  elevator 
input. 

The  last  but  most  important  improvement  has  already  been  dis- 
cussed in  some  detail;  the  bomb  should  be  moved  further  from  the  air- 
craft center  of  gravity.     The  importance  of  this  change  is  shown  in 
Figs.  30  through  33  which  present  analog  computer  time  histories  of 
pitch  rate  and  normal  acceleration.     The  moment  input  here  corresponds 
to  a  distance  of   1.0  feet  between  the  aircraft  and  bomb  centers  of  gravity, 
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The  effect  of  each  stability  derivative  on  thi.6  type  of  response  is  shown. 
In  these  comparisons,  the  same  numerical  values  were  used  as  those 
employed  with  the  flight  test  comparisons  shown  in  Figs.  20  through  24. 
With  larger  distance  from  the  bomb  to  the  aircraft  e.g.  ,  the  moment 
input  could  be  determined  accurately.     These  curves  show  that  each 
stability  derivative  has  a  larger  effect  on  the  transient  response  when  a 
larger  moment  input  is  used.     The  relative  importance  of  each  deriva- 
tive  is  not  greatly  effected; -however        "Ma;      appears  to  show  the 
greatest  improvement. 


There  is  a  characteristic  of  these  curves  which  is  of  consider- 
able assistance  in  determining  the  derivatives  by  curve  matching. 
There  is  a  short  time  where        ^a      and        "a      have  no  effect  on  the 

transient  response.     This  is  represented  on  the  curves  by  the  crossing 

C  CT 

of  the  variables    n    and     dO     as        mo.'      and        J-'  a     are  varied. 

This  means  physically  that  the  aircraft  is  passing  through  zero  angle  of 

attack.     Fig.   34  shows  that  angle  of  attack  initially  decreases  and  then 

increases  through  zero  to  a  positive  steady  state  value.     In  the  vicinity 

C 

of  the  minimum.        mda      has  no  effect  on  the  transient  response  while 

c  c 

m^      and        La      have  no  effect  where     a      passes  through  zero.     In 

the  flight  tests,  the  steady  state  angle  of  attack  was  negative  and  no 
zero  value  occurred. 


If  the  forcing  functions  are  known  accurately,  approximate 
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values  for  the  stability  derivatives  will  locate  the  time  periods  where 

0!     is  very  small  and  where  the  minimum  occurs    {da    -    0).      These 

C  C 

regions  may  then  be  used  to  determine        radQ     and        mda      independ 

C  C 

ent  of         La     and       ma  . 


Ct 

The  use  of  an  elevator  step  to  determine        ^a      initially  is 

advantageous.    Since  an  angle  of  attack  vane  can  be  used  with  this  input, 

CT 
time  histories  of     n     and     a      give        ■Ljq!      directly.     For  the  elevator 

3tep,  equations  (7)  and  (13)   show  that: 


V  2    T 


a 


With        ^a      known,  the  remaining  derivatives  can  be  determined  with 
the  analog  computer  as  follows.    If  a  reasonably  good  match  is  not  avail  - 

CT 

able  initially,  vary  all  stability  derivatives  except        ^a      to  achieve  a 
fair  correlation  throughout.      Leave        mda      fixed  and  vary         md9 

and        ma      to  match  the  regions  where     da      is  very  small  (the  mini- 

cm 
mum  and  the  steady  state  of    a  ).      Then  vary  da      to  match  the 

initial  portions  of  the  time  histories. 


When  a  good  match  is  achieved,   shift  to  the  elevator  step  input. 

c  c  c  c 

Leave        mdQ!      and        md0       fixed  and  adjust        m  a      and        La      to 
achieve  the  best  match  for  all  three  variables    (n,  a    and    dG).    Return 
to  the  bomb  input  and  leave        ^a     and        IJla      fixed  and  achieve  the 
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C  C 

best  match  possible  varying        mda      and        mdQ  .      Shift  between  the 

two  forcing  functions  using  these  techniques  until  the  best  over-all 

match  is  achieved. 

While  the  bomb  drop  establishes  all  derivatives  uniquely,  the 

use  of  the  elevator  step  will  considerably  improve  the  accuracy  with 

C—j  Ct 

which  a      and  a      can  be  determined.     The  accuracy  of  the 

rate  derivatives  is  improved  only  slightly.  Since  an  angle  of  attack 
vane  can  be  used  successfully  with  an  elevator  step,  an  extra  cross 
check  is  available. 
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CONCLUSIONS 

The  proposed  method  of  uniquely  determining  longitudinal 
stability  derivatives  for  slow  speed  aircraft  is  believed  to  be  feasible. 
Due  to  instrumentation  difficulties,   positive  verification  of  this  is  lack- 
ing.    However,   the  numerical  results  for  the  Navion  contained  in  Table 
III  are  considered  reasonably  valid. 

The  stability  derivatives  were  determined  uniquely  from  simu- 
lated data  generated  on  the  analog  computer .     The  simplest  and  most 
accurate  results  were  obtained  by  a  curve  matching  technique. 

The  use  of  an  elevator  step  function  in  addition  to  a  weight 
release  will  greatly  assist  in  the  determination  of  the  longitudinal  sta- 
bility derivatives.     No  consideration  was  given  to  any  other  types   of 
forcing  functions. 

The  accuracy  with  which  each  stability  derivative  is  determined 
is  associated  with  its  influence  on  the  aircraft  response.     The  derivative 

Q 

most  accurately  determined  by  the  release  of  a  weight  is  x   du  > 

C_ 

da      is  uniquely  found  with  fair  accuracy  by  a  bomb   release. 

GL      and         ^ct      can  be  found  with  better  accuracy  using  an  elevator 

m  m 

step  and  previously  determined  values  for  dO       and  dc  . 
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RECOMMENDATIONS 

As  a  result  of  this  investigation,  the  following  recommenda- 
tions are  proposed: 

1)  An  investigation  of  this  nature  should  be  continued  to  obtain 
improved  experimental  results.  The  configuration  and  the  velocity  of 
the  aircraft  should  be  varied  in  any  additional  investigation. 

2)  A  normal  accelerorneter  with  a  full  scale  range  of  approxi- 
mately one     g     should  be  installed  in  an  aircraft  for  longitudinal 
transient  flight  testing. 

3)  The  aircraft  should  be  instrumented  to  measure  elevator 
position  and  attitude  to  properly  analyse  the  response  from  a  bomb 
release. 

4)  In  this  type  of  investigation  in  a  Navion  aircraft,  a  500  pound 
weight  should  be  at  least  one  foot  horizontally  from  the  aircraft  center 
of  gravity. 
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TABLE  I 


Physical   Characteristics   of  the   Navion 
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A.      Wing 


1.  Total  Wing  Area  (includes  flaps,   ailerons 

and    19.8  ft2   of  fuselage) 

2 .  Span 

3.  Mean  Aerodynamic  Chord 


4. 

Angle  of  Incidence 
Root 
Tip 

5. 

Twist 

Aerodynamic 
Geometric 

6. 

Airfoil  Section 
Root 
Tip 

7. 

Aspect  Ratio 

8. 

Taper   Ratio 

9. 

Dihedral 

10. 

Chord 

Root 
Tip 

184.2  ft2 

3  3.38  ft 

5.70  ft 

+2° 

-10 

2°  31' 

3° 

NACA  4415R 

NACA  641  OR 

6.044 

0.5265 

7.5° 

7.2   ft 

3 .  92  ft 

Vertical   Tail 


1.      Total  Area 


2.      Stabilizer   Area 


3.      Rudder  Area 


12.93  ft2 
6.87  ft2 
6.05  ft2 
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TABLE   I    (Cont.) 


4.  Airfoil  Section 

Root 

Tip 

5.  Angle  of  Stabilizer  with  fuselage 


NACA  0013.2     Mod, 
NACA  0012-64  Mod, 

2°  Nose  left 


C.     Horizontal   Tail 

1.  Total  Area 

2.  Stabilizer   Area  , 

3.  Elevator   Area 

4 .  Span 

5.  Mean  Aerodynamic  Chord 

6.  Angle  of  Incidence 

7.  Airfoil  Section 

8.  Aspect  Ratio 

9.  Taper   Ratio 


43.05     ft 
28.95     ft2 
14.10     ft2 

13.  172  ft2 
3 .  34     ft2 

-4° 

NACA  0012 

4.02 

.67 


D. 

General 

1. 

Length  Overall 

2. 

Tail   Length 

3. 

Weight    (full  fuel) 

4. 

Engine 

27.25  ft 

16.88  ft 

2430  pounds 

Continental  E-185 
185  HP  at  2300  RPM, 
29"  HgM.P.  at  sea  level 


TABLE  II 


Test   Flight  Conditions 
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Flight  4 

Flight  9 

Flight  9B 

Date 

23  April  59 

2  May  59 

2  May  59 

Airspeed  mph 

94 

.     96 

94 

Temp 

7°C. 

9°C. 

9°C. 

Pressure 

Altitude   (ft) 

350 

150 

450 

e.g. 

29 . 2  %  mac 

32.  5%  mac 

32.  5%  mac 

Vi  (fps) 

137.8 

140.5 

137.8 

V    (fps) 

136.6 

139.4 

137.3 

or 

1.018 

1.0175 

1.006 

w 

2540 

2538 

2525 

CL 

.748 

.720 

- 

CL   ' 
^O 

.611 

.589 

.607 

ACL 

.  1360 

.1312 

I 

3132 

3541 

3541 

h 

.0787 

.0887 

.0887 

T 

1.296 

1.273 

1.298 

TABLE  in 


Stability  Derivatives  from  Flight  Test6 
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Flight  4 

Flight  9  and  9B 

Date 

23  April  59 

2  May  59 

e.g.  locati 

on 

29.  2%  mac 

32 .  5%  mac 

CL 

4„98 

4.98 

cma 

-.325- 

-  .60 

CmdG! 

• 

-  .072 

-  .073 

Cmd0 

-  .167 

-  .166 

vm 

.00244 

.00463  (Flight  9,  bomb  drop) 

Calculated 

A 

Cm 

.0029 

.0065 

for   Bomb  Drop 
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Fig.   2       Bomb  Rack  Installation 


Fig.   3        Test   Aircraft 


Fig.   4        Tape  Recorder  and  Telemeter  Receiver 


Fig.    5       Analog  Computer  and  Sanborn  Recorder 
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